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Abstract 

This comment is in response to the paper by L. Guzzo recently ap- 
peared in "New Astronomy" related to our work. The subject of the 
discussion concerns the correlation properties of galaxy distribution in 
the available 3-d samples. There is a general agreement that galaxy 
structures exhibit fractal properties, at least up to some scale. How- 
ever the presence of an eventual crossover towards homogenization, as 
well as the exact value of the fractal dimension, are still matter of de- 
bate. Here we briefly summarize our point of view by discussing three 
main topics. The first one is methodological, i.e. we clarify which are 
the correct methods to detect the real correlations properties of the 3-d 
galaxy distribution. Then we discuss the results of the analysis of sev- 
eral samples in two ranges of scales. In the first range of scale, below 
100 -v- 200/i~ 1 Mpc, the statistical quality of the data is rather good, and 
we find that galaxy distribution has fractal properties with D ~ 2. At 
larger distances the statistical robustness of the present data is weaker, 
but, we show that there is evidence for a continuation of the fractal 
behavior, without any tendency towards homogenization. 



1 Introduction 



In a recent paper L. Guzzo |§ has exposed his arguments in favor of the 
homogeneity of galaxy distribution in the available three dimensional sam- 
ples. This paper takes its origin from a discussion between Dr. Guzzo and 
one of us (F.S.L.), held during the fifth Italian National Cosmology Meet- 
ing (Dec. 1996). This debate followed the one occurred in the Conference 
"Critical Dialogues in Cosmology" between Prof. M. Davis and L.P. (June 



1996) (see |2]J and |4[). Moreover as a supplement to the debate between 
L.P. and M. Davis, Prof. Peebles has sent a circular letter where he exposes 
his arguments in favor of homogeneity, followed by a similar letter that L.P. 
sent to Prof. Peebles. All this material is now available at the home page: 
http://www.phys.uniromal . it/DUCS/PIL /pit. htm\ 

In this comment we briefly summarize, in a colloquial way, our opinion 
about the arguments of Dr. Guzzo and we refer the reader to a recent review 



28] for a more comprehensive and detailed discussion of the subject. 

Dr. Guzzo puts a special emphasis on the fact that the identification of the 
scale at which galaxy distribution becomes homogeneous is one of the major 
topics of modern Cosmology. Homogeneity is, in fact, the basic assumption of 
the current theory of galaxy formation, and in general is the cornerstone of any 
cosmological theory. Moreover it has been elevated to the role of a principle by 
the Cosmological Principle (e.g. We have discussed in several papers || 

23 1 p8fl that homogeneity is a very strong a-priori assumption, although it has 



been a quite reasonable one up to the compilation of large redshift surveys. As 
Mandelbrot stressed several times |lTJ [|1^], the assumption of local isotropy, 



without analyticity of matter distribution, ensures as well the equivalence of 
all the occupied points in the universe. Hence from a conceptual point of view, 
an isotropic fractal distribution is completely compatible with the requirement 
of having no special directions or positions in the universe. 

Therefore, it is a crucial challenge in nowadays Cosmology to test whether 
the matter distribution is analytical or not, or in other words, whether there 
exists a characteristic length in galaxy (and cluster) distributions. To this end 
a particular care must be used in the discussion of the correlation properties 
of galaxy samples. As Guzzo does, we focus here on the properties of redshift 



samples, and we refer the reader to [28] [14 and for a more detailed discussions 



of the angular properties of real distributions, and to [JTT1 [0] (T^[ JTJj for those 
of artificial fractals. 

Dr. Guzzo correctly pointed out that "the specific scale at which the galaxy 
distribution apparently turns to homogeneity is dangerously close to the size 
of the largest sample presently available", i.e. ~ 100 -j- 200h~ 1 Mpc. However 
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he claims that at larger distances there are enough evidences to conclude that 
homogeneity is well established. Here we will argue that there are well defined 
and ample evidences that the transition scale to homogeneity is larger than 
100-=-200/i _1 Mpc, and that up to this distance galaxy (and cluster) distribution 
shows well defined power law correlations, corresponding to a fractal dimension 
D sa 2. Hence we partially agree with Dr. Guzzo about the small scale 
properties of galaxy distribution. However, we find in his arguments a lack 
of discussion of the implications of this result, especially for what concerns 
the luminosity segregation effect and, more in general, for the methods used 
to characterize galaxy correlations at these distances. For example if galaxy 
distribution becomes homogeneous at ~ 100h~ 1 Mpc, what is the meaning of 
the "correlation lengths" of galaxies and clusters ? These are important points 
having a number of implication, which we consider here in more detail. 

Moreover we discuss our point of view, about galaxy and cluster distribu- 
tion at scales larger than 100 4- 200h~ 1 Mpc, trying to clarify the statistical 
robustness of our results. In this range of length scales we disagree with Dr. 
Guzzo. In particular we show that, although the statistical quality of the data 
at these distances is weaker than at smaller scales, there is no any evidence 
for homogenization in any of the 3-d samples published up to now, and on 
the contrary there are evidence which support the continuation of the fractal 
behavior with dimension D « 2 found at smaller distances. 



2 A methodological point 

There is a general agreement that galaxy distribution exhibits fractal behavior 



up to a certain scale (e.g. ]TT[ [17j [[J |18fl ). The eventual presence of a transi- 



tion scale towards homogenization and the exact value of the fractal dimension 



are matters of the present debate |2T| Q . Given this situation, it is first of all 
crucial to establish which are the statistical methods suitable and appropriate 
to characterize the correlation properties of galaxy distribution. In particu- 
lar we answer to the following question: which are the statistical tools able 
to eventually identify the homogeneity scale and to measure the correlation 
exponent in the correct way ? 

The proper methods to characterize irregular as well as regular distributions 
have been correctly illustrated by Dr. Guzzo in Sec. 3, and we refer to |§ and 



28] for a more detailed and exhaustive discussion. The basic point is that, 
as far as a system shows power law correlations, the usual £(r) analysis (e.g. 
P7| ) gives an incorrect result, since it is based on the a-priori assumption of 



homogeneity. In order to check whether homogeneity is present in a given 
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sample one has to use the conditional density T(r) defined as JHJ 

, , (n(rS)n(r* + r)) AD n , 
(n) 47r 

where the last equality holds in the case of a fractal distribution with dimension 
D and prefactor A (we follow here the same definitions of ||). In the case of an 
homogenous distribution (D = 3) the conditional density equals the average 
density in the sample. Hence the conditional density is the suitable statistical 
tool to identify fractal properties (i.e. power law correlations with codimension 
7 = 3-D) as well as homogeneous ones (constant density with sample size). If 
there exists a transition scale Ao towards homogenization, we should find T(r) 
constant for scales r > A . 

It is simple to show that in the case of a fractal distribution the usual £(r) 
function in a spherical sample of radius R s is [|(J || 



D ( r \ D ~ 3 

From Eq.||] we can see two main problems of the £(r) function: its amplitude 
depends on the sample size R s (and the so-called correlation length r , defined 
as £( r o) = 1, linearly depends on R s ) and it has not a power law behavior. 
Rather the power law behavior is present only at scales r < ro, and then it 
is followed by a sharp break in the log-log plot as soon as £(r) ^ 1. Such a 
behavior does not correspond to any real change of the correlation properties 
of the system (that is scale - invariant by definition) and it makes extremely 
difficult the estimation of the correct fractal dimension as it shown in Fig.|l|. In 
particular if the sample size is not large enough with respect to the actual value 
of r , the codimension estimated by the £(r) function (7 w 1.7) is systematically 
larger than 3 - D (7 re 1) |28fl. 

Given this situation it is clear that the £(r) analysis is not suitable to be 
applied unless a clear cut-off towards homogenization is present in the samples 
analyzed. As this is not the also Dr. Guzzo stressed, it is appropriate 



and convenient to use T(r) instead of £(r). We have discussed in detail in |28 
that the use of the correct statistical methods is complementary to a change 
of perspective from a theoretical point of view. 



3 Galaxy distribution at scale r < 100^-200/i _1 M^ c 

From the previous discussion it seems that Dr. Guzzo agrees with us about the 
fact that, unless a well defined value of the average density has been established, 
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Figure 1: The £(r) function computed for a fractal with dimension D = 2 
(7 = 3 — D = 1) in two samples with different depth R s . In the first case 
R s = 20Mpc one obtains an higher value for the correlation exponent 7 = 1.7 
because one performs the fit in the region of length scales r ~ ro, as the power 
law behavior is poorly defined. In the second case the power law behavior 
of £(r) is more extended (R s = lOOMpc) and it is now possible to fit in the 
region r<Cr . Here one obtains the correct value for the correlation exponent 
7 = 1. However the effective depth of the various galaxy samples is in the 
range ~ 10 -7- 30/i _1 Mpc (see Table 1) and hence this effect is very important 
in the real data. 
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Table 1: The volume limited catalogues are characterized by the following 
parameters: - Rai^h' 1 Mpc) is the depth of the catalogue - Q is the solid angle 
- R s (h~ 1 Mpc) is the radius of the largest sphere that can be contained in the 
catalogue volume. This gives the limit of statistical validity of the sample. - 
ro(h~ 1 Mpc) is the length at which £(r) = 1. - Aq is the eventual real crossover 
to a homogeneous distribution that is actually never observed. The value of 
r is the one obtained in the deepest VL sample. The CfA2 and SSRS2 data 
are not yet available and in this case the values of r have been taken by. 
Park et al. 1994 and Benoist et al. 1996 respectively. (Distance are expressed in 
h^Mpc). 



Sample 


fi (sr) 


Rd 


R s 




D 


Ac 


CfAl 


1.83 


80 


20 


6 


1.7 ±0.2 


> 80 


CfA2 


1.23 


130 


30 


10 


~ 2.0 


? 


PP 


0.9 


130 


30 


10 


2.0 ±0.1 


> 130 


SSRS1 


1.75 


120 


35 


12 


2.0 ±0.1 


> 120 


SSRS2 


1.13 


150 


50 


15 


~ 2.0 


? 


Stromlo-APM 


1.3 


100 


35 


12 


2.2 ±0.1 


> 150 


LEDA 


~ 5 


300 


150 


45 


2.1 ±0.2 


> 150 


LCRS 


0.12 


500 


18 


6 


1.8 ±0.2 


> 500 


IRAS2 Jy 


~ 5 


60 


20 


5 


2.0 ±0.1 


> 50 


IRAS1.2Jj/ 


~ 5 


80 


30 


8 


2.0 ±0.1 


> 50 


ESP 


0.006 


700 


8 


3 


1.9 ±0.2 


> 700 



and found to be independent on the sample size, the usual £(r) analysis fails. 
Hence we should focus on the determination of the average density rather than 
on r : the latter quantity being meaningful unless an homogeneous distribution 
has been found. 

The first test on fractal versus homogeneous properties concerns the relation 
between the sample size R s and the so-called correlation length tq measured in 
redshift catalogs. The conclusion of Dr. Guzzo is that this relation is not the 
linear one predicted in the case of fractal distributions. Here we briefly report 
our analysis of the tq — R s relation, which significantly disagrees with the one 
of Dr. Guzzo. 

In Tab.|l| we report the characteristics of the various catalogs we have an- 
alyzed by using the methods illustrated in Sec. 2. Here we consider in more 
detail only three catalogs, where our estimations of r and/or R s significantly 
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disagree with those of Guzzo and a detailed explanation of the analyses of the 
catalogs shown in Tab.l can be found in [pgfl . However before doing this we 
would like to stress two important points: 

1. Given a certain sample of solid angle Q and depth Rd, it is important 
to define which is the maximum distance up to which it is possible to compute 
the correlation function (T(r) or £(r)). As discussed in |J, we have limited 
our analysis to an effective depth R s that is of the order of the radius of 
the maximum sphere fully contained in the sample volume. In such a way 
we eliminate from the statistics the points for which a sphere of radius r is 
not fully included within the sample boundaries. Hence we do not make any 
assumption on the treatment of the boundaries conditions. Of course, by doing 
this, we have a smaller number of points and we stop our analysis at a smaller 
depth than that of other authors. 

The reason why T(r) (or £(r)) cannot be computed for r > R s is essen- 
tially the following. When one evaluates the correlation function (or the power 
spectrum [B7[]) beyond R s , then one makes explicit assumptions on what lies 
beyond the sample's boundary. In fact, even in absence of corrections for se- 
lection effects, one is forced to consider incomplete shells calculating T(r) for 
r > R s , thereby implicitly assuming that what it is not seen in the part of 
the shell not included in the sample is equal to what is inside (or other simi- 
lar weighting schemes). In other words, the standard calculation introduces a 
spurious homogenization which we are trying to remove || [ 28|| . 



We have done a test [EHf on the homogenization effects of the incomplete 



shells on artificial distributions as well as on real catalogs, finding that the 
flattening of the conditional density is indeed introduced owing to the weight- 
ing, and does not correspond to any real feature in the galaxy distribution. 
These results differ from those of p2| (this is reported in the Appendix A in 
Dr. Guzzo's paper) probably because they did not take into account finite 
size effects in the generation of artificial samples and they considered ensemble 
averages of the conditional density (see for a more detailed discussion). 

2. We do not use weighting schemes, and hence our analysis concerns only 
volume limited (VL) samples. The use of magnitude limited (ML) samples 
and the weighting schemes inevitably requires a-priori assumptions on nature 
of the distribution |28| . 



Now we discuss in detail our disagreement with Tab.l of the paper of Guzzo. 

• ESP. In this case the estimation of R s slightly differs from that of Guzzo, 
probably because we have not used the relativistic corrections ( [PS|1 see 
below). Also the value of ro is slightly different (tq ~ 3/i~ 1 Mpc instead of 
the measured r ~ 4.5/i _1 Mpc), and this is probably due to the fact that 
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ESP does not cover a continue solid angle in the sky, as it is a collection 
of pencil beams. Such a situation necessarily requires the introduction of 
spurious treatments of the boundary conditions (see point 1). 

LCRS. This survey has the peculiar property of being limited by two 
limits in apparent magnitude (a lower and an upper one). In order to 
construct a VL sample in this case, one has to impose two limits in 
distances and correspondingly two in absolute magnitude. This is the 
origin of a smaller R s in our Table 1 than this reported by Guzzo (R s ~ 
32/i _1 Mpc) ||28|| . This implies a smaller r , much closer to the measured 
one. 

Stromlo/APM. We have extensively analyzed this catalog in and 
|29| and the value of r is reported in Tab.l. Due to the sparse sampling 
strategy adopted to construct this catalog, we are able to measure the 
correlation properties up to R s ~ 40h~ l Mpc and not 83h~ 1 Mpc as re- 



ported by Guzzo. The disagreement with the work of Loveday et al. jlO 
(ro ~ Y2h~ x Mpc rather than ro « 5/i _1 Mpc) is probably due to the treat- 
ment of the boundary conditions and their use of ML samples rather than 
VL ones (i.e. they used weighting schemes with the luminosity selection 
function). In any case, we stress again, the proper test is check whether 
the conditional density has a power behavior [2~9~|. 



We show in Fig.^| the results of the conditional density determinations in 
various redshift surveys |23| [E5|. All the available data are consistent with 
each other and show fractal correlations with dimension D = 2.0 ± 0.2 up 
to the deepest scale probed up to now by the available redshift surveys, i.e. 
~ 150ft -1 Mpc. A similar result has been obtained by the analysis of galaxy 



cluster catalogs [|T§] [p8 



Finally it must be noted that Dr. Guzzo based a part of his arguments 
on the fact that a luminosity bias is responsible of the shift of ro with sample 
size (see also [0 We do not enter into a detailed explanation of the 

inconsistencies of the usual argument about luminosity segregation |L6| 



and we refer to [28| for a more detailed analysis. However we would like to 



remark that the authors (e.g. [16] 0) who have addressed this concept have 
never presented any quantitative argument that explains the shift of r with 
sample size. In this respect, an exception is represented by the paper of Davis 
et al. (1988) in which the authors claim that ro behaves as R° s 5 . However for 
what concerns luminosity segregation, the meaningful parameter must be the 
absolute magnitude limit of the volume limited (VL) sample considered rather 
than the depth (R s ). Having fixed the limiting apparent magnitude of the 
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Figure 2: The spatial density T(r) computed in some VL samples of CfAl, 
PP, LEDA, APM, ESP, LCRS, SSRS1, IRAS and ESP (Sylos Labini et al, 
1997). 
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catalog, at each R s would correspond a well defined absolute magnitude. The 
brightest galaxies are present yet in samples like CfAl, and hence according 
to the luminosity segregation paradigm, there is no reason one should expect 
that in deeper sample (like CfA2 or SSRS2) r$ is increased. However this is 
actually the case [ |T6|j . 



We have discussed in detail in [24] that the observation that the giant 



galaxies are more clustered than the dwarf ones, i.e. that the massive elliptical 
galaxies lie in the peaks of the density field, is a consequence of the self-similar 
behavior of the whole matter distribution. The increasing of the correlation 
length of the £(r) has nothing to do with this effect, rather it is related to the 
sample size. 

A last point requires a further clarification. Dr. Guzzo claims that T(r) is a 
power law with D ~ 1.2 up to r ~ 3.5h~ 1 Mpc, then it shows a different scaling 
between ~ "ih^Mpc and 30h~ 1 Mpc with fractal dimension D ~ 2 ||. In our 
view there is a subtle point in this discussion that has not been previously 
appreciated. 

Suppose now, for simplicity, we have a spherical sample of volume V in 
which there are iV points, and we want to measure the conditional density. It 
is possible to compute the average distance between neighbor galaxies (A), in 
a fractal distribution with dimension D, and the result is 

(A)= (\Yt(i + ^) (3) 



A) \ D 

where V is the Euler's gamma-function (Note that the prefactor A is 

dependent on the luminosity selection function of the VL chosen). Clearly 
this quantity is related to the lower cut-off of the distribution A (eq.|IJ) and to 
the fractal dimension D. If we measure the conditional density at distances 
r < (A), we are affected by a finite size effect. In fact, due the depletion of 
points at these distances we underestimate the real conditional density finding 
an higher value for the correlation exponent (and hence a lower value for the 
fractal dimension). In the limiting case at the distances r <C (A), we can 
find almost no points and the slope is 7 = —3 (D = 0). In general, when 
one measures T(r) at distances which correspond to a fraction of (A), one 
finds systematically an higher value of the conditional density exponent. Such 
a trend is completely spurious and due to the depletion of points at such 
distances. It is worth to notice that this effect gives rise to a curved behavior of 
T*(r) (the integral of T(r) see || fl2"8| ) at small distances, because of its integral 



nature. This is exactly the case of the deepest VL of Perseus-Pisces which 
Guzzo et al. (1991) considered in their analysis, and for which (A) ~ 8/i _1 Mpc. 
A clarifying test in this respect would be to check whether this change of slope 



10 



is actually present also in the others VL samples of the same survey, which 
have a larger number of points (and hence a lower (A)). This test has been 
performed by [[25] and the conclusion is that the change of slope is due a finite 
size effect rather being an intrinsic property of galaxy distribution. 

Finally we would like to point out another point that is inconsistent in the 
arguments of Dr. Guzzo. He and his collaborators found that T(r) flattens 
at > 30/i~ 1 Mpc. Up to this distance r must be a linear fraction of R s in view of 
eq|J and does not depend to any luminosity bias ! In our opinion [^6] J^J , this 
flattening is due to an incorrect treatment of the boundary conditions (see point 
1). However if this behavior is real, Dr. Guzzo should conclude that galaxy 
distribution is homogeneous at ~ 30h~ 1 Mpc and not at ~ 100 -j- 200/zT 1 Mpc 
as he claims. 

We conclude this discussion pointing out another element: the analysis of 
LEDA. Dr. Guzzo claims that the results coming from LEDA B fTTTl 
having no meaning whatsoever" . In our opinion such a strong conclusion must 
be supported by quantitative arguments. We have done several tests on the 
LEDA database || IJ and we have concluded that, although this sample is 
highly incomplete, the statistical results obtained are rather stable and robust. 
We refer to the previously mentioned papers for a more specific discussion. 



4 Galaxy distribution at scale r > 100^-200/i 1 Mpc 

We discuss now the behavior of the radial density in general, and then we 
consider the two cases of PP and ESP (see [T3] and [2S| for a more detailed 
discussion on this point). 

We focus on the possibility of extending the sample effective depth R s . In 
order to discuss this question, it is important to analyze the properties of the 
small scale fluctuations. To this aim, we introduce the conditional density in 
the volume V(r) (that can be a portion of a sphere) as observed from the origin, 
defined as 

^^^-3. (4 ) 

where the factor p comes from the fact that a VL sample contains only a fraction 
Ny L (R) = p ■ N(< R) (where p < 1) of the total number N(< R) of galaxies 
in V(R). If 4>(L)dL is the fraction of galaxies whose absolute luminosity (L) is 
between L and L + dL, p is given: 

JV 00 <P(L)dL 

'■' 'nil) 1 v J 
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In Eq.[| Lvl is the minimal absolute luminosity that characterizes the VL sam- 
ple and L min is the fainter absolute luminosity (or magnitude M min ) surveyed 
in the catalog (usually M min ~ —11). Computing Ti(r), we expect (Fig.|3]) not 
to see any galaxy up to a certain distance (A). For distances somewhat larger 
than (A) , we expect therefore a raise of the conditional density because we are 
beginning to count some galaxies and Tj(r) is affected by the fluctuations due 
to the low statistics. 

It is therefore important to be able to estimate and control the minimal 
statistical length A, which separates the fluctuations due to the low statistics 
from the genuine behavior of the distribution. A simple argument for the 
determination on the length A is the following [25[] ||28|| . At small scale, 



where there is a small number of galaxies, there is an additional term, due 
to shot noise, superimposed to the power law behavior of Tj(r), that destroys 
the genuine correlations of the system. ( As we have discussed in |2S| there 



should be considered also an intrinsic oscillating term for non-averaged out 
quantities: however for sake of clarity we avoid this discussion here.) Such a 
fluctuating term can be erased out by making an average over all the points 
in the survey. On the contrary, in the observation from one point, when the 
number of galaxies is large enough the shot noise becomes negligible. Roughly, 
this happens when the number of points is lager than ~ 30 (see Fig|J). This 
condition gives (from Eq.|j) 

/ 4tt \* 20 -r- fflhr l Mvc 

for a typical VL sample with Myi ~ M*, where A corresponds to the amplitude 
of the conditional density of all galaxies [25| f28|. This can be estimated from 



the amplitude of T(r) in a VL sample divided by the correspondent p as defined 
in Eq.|5]. We find (for typical catalogues) B « 10 -f- 15(/i" 1 Mpc) _ ' D f25| . The 
corresponding value of A are reported in Tab. 2. 

The difference between T(r) and r,(r) is straightforward: the latter quan- 
tity is not an average one. However it is an integrated one. The differential 
density is shells is clearly much more noisy, because it is not averaged neither 
integrated. We have computed the Tj(r) in the various VL samples of Perseus- 
Pisces redshift survey, and we show the results in Fig.[|. In the less deeper VL 
samples (VL70, VL90) the smooth behavior of the density is weakly defined 
because in this case the finite size effects are very important as the distances in- 
volved are r < A (Eq.|]). We note that about the same scales we find a very well 
defined power law behavior by the correlation function analysis. In the deeper 
VL samples (VL110, VL130) a smooth and well defined behavior is reached 
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Figure 3: Behavior of the density computed from one point, in the case of 
a fractal distribution. At small distances below the average mean separation 
between neighbor galaxies, one finds no galaxies. Then the number of galaxies 
starts to grow, but this regime is strongly affected by finite size fluctuations. 
Finally the correct scaling region r m A is reached. In the intermediate region 
the density can be approximated roughly by a constant value. This leads to 
an apparent exponent D ps 3 for the integrated counts. This exponent is not 
real but just due to the size effects. 
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SllTA/'PA/' 


Of sr) 


Affr _1 /\/fwr) 


CfAl 


1.8 


15 


CfA2 (North) 


1.3 


20 


SSRS1 


1.75 


15 


SSRS2 


1.13 


20 


PP 


1 


40 


LEDA 


2 7T 


10 


IRAS 


2 7T 


15 


ESP 


0.006 


300 



Table 2: The minimal statistical length A for several redshift surveys. 



for distances larger than the scaling distance (Q = 0.9 sr) r w A ~ 50/i" 1 Mpc. 
The fractal dimension is D ~ 2 as the one measured by the average conditional 
density. 

For relatively small volumes it is possible to recover the correct scaling 
behavior for scales of order of (A) (instead of ~ 10(A)) by averaging over several 
samples or, as it happens in real cases, over several points of the same sample 
when this is possible. Indeed, when we compute the correlation function, we 
perform an average over all the points of the system even if the VL sample is 
not deep enough to satisfy the condition expressed by Eq.|. 

Now, if we look at the density in shells, this will be dominated by fluc- 
tuations, as it is not a cumulative distribution. Note that if the distribution 
would became homogenous, as Guzzo et a/. found, at ~ 30/i _1 Mpc, then at 
3 times the homogeneity scale the distribution would present a rather regular 
behavior. This is clearly not the case. 

Let us now briefly discuss the properties of the ESP catalog. As the redshifts 
involved are quite large a particular care is devoted to the construction of the 
VL subsamples for the case of ESP. We have used various distance-redshift 
(d(z)) and magnitude-redshift (m(z)) relations. Moreover the data are selected 
in the blue-green and even if the redshifts of galaxies are moderate (z < 0.2), 
K-corrections are needed to compute the absolute magnitude of galaxies. The 
corresponding magnitude-redshift relation is 

m-M = 25 + 5 log 10 (d(z, q )(l + z)) + a(T)z (7) 

where we have used the functional forms of the K-correction a(T)z as a function 
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PERSEUS-PISCES 



• PP-130 




» PP-70 




10 100 
r (Mpc) 

Figure 4: The spatial density Tj(r) computed in the VL sample cut at 
70,90, 110, UOh^Mpc . In the case of VL70 and VL90 the density is domi- 
nated by large fluctuations and it has not reached the scaling regime. In the 
samples VL110 and VL130 the density is dominated by large fluctuations only 
at small distances, while at larger distances, after the Perseus Pisces chain at 
a very well defined power law behavior is shown, with the same 
exponent of the correlation function (i.e. D = 2) 
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of redshift. a(T) depends on the morphological type T and goes from a ~ 2 for 
the Scd galaxies to a ~ 3.7 for the E/SO galaxies. It is not possible to apply the 
K-correction to each morphological type because over the VI th magnitude it is 
not possible to recognize the Hubble type from visual inspection. To overcame 
this problem we have adopted a statistical approach: we have assumed various 
percentage of late and early-type galaxies. The observed percentage in nearer 
samples is 70% of late and 30% of early type galaxies. We have performed 
a number of tests by varying these percentages to show that the final result 
depends weakly on the adopted values. We note that varying these percentages 
change the number of points in the VL subsamples as the absolute magnitude 
can change of about a unit. 

The behavior of the number- distance relation in the standard FRW model 
depends strongly on the value of the deceleration parameter q Q for high redshift 
z > 0.5, while for z < 0.2 the relativistic corrections are very small for all 
reasonable q . 

All the samples show a highly fluctuating behavior for the space density as 
it is reported by Dr. Guzzo. It is very difficult in this case to identify a clear 
power law behavior. However we stress that an homogeneous distribution (say 
at 50h~ 1 Mpc) would show a very smooth and flat behavior for Ti(r) at these 
distance scales. This is clearly not the case. 

It is very interesting, in our opinion, to study also the case of a purely 
Euclidean space. In this case we can write the d(z) relation simply as: 



and the m(z) relation (without K-corrections): 

m-M = 5 log 10 + 25 . (9) 

In this case the volume grows as ~ z 3 ~ d 3 . The behavior of the radial density 
in the Euclidean case is reported in Fig.[|[ From this figure it is possible 
to see the effect of the finite size fluctuations: the scaling region begins at 
A ~ 300 -T- 320h^ 1 Mpc for all the VL subsamples. For smaller distances the 
statistical fluctuations dominate the behavior of Tj(r). It seems that in this 
case the power law behavior for r > A ~ 320h~ 1 Mpc is better defined than in 
the previous case. The fractal dimension turns out to be D w 2.2. 

The K-correction, that must be applied for the Hubble shift of the galaxy 
spectra, can change the absolute magnitude of a unit. This correction is due 
to a systematic effect for each morphological type. As we put these correction 
in a random way we are mixing a systematic effect with a random correction. 
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Figure 5: The spatial density Tj(r) computed in a VL samples of ESP. A the 
power law behavior with D pa 2 is shown for r > 300h~~ 1 Mpc ~ A, even though 
with strong fluctuations 



17 



In this way we can only check the stability of the results but we cannot hope 
to obtain a better fit p8f By doing this we find a marginal power law 
decaying of the conditional density. Therefore we may conclude that there is 
a weak evidence that the fractal dimension is D « 2 in this sample due to the 
poor statistics. 

There are several other evidences in our opinion that point towards a fractal 
distribution of galaxies at very large scale, and in particular they are the density 
behavior from one point in the LCRS and the behavior of the galaxy numbers 
counts as a function of the apparent magnitude [p8j| . However we stress again 
that, due to the lackness of complete redshift measurements, these evidences 
are statistically weaker than the ones up to one hundred Megaparsec. 



5 Discussion and Conclusions 

Only after clarification of the small scale galaxy correlations it is possible to 
investigate the large scale ones. For what concerns the meaning of the so-called 
"correlation lengths" of galaxies and clusters, the behavior of the conditional 
density up to ~ 100h~ 1 Mpc is enough to give us the elements for a revision 
of both the statistical methods usually used in the data analysis as well as 
the theoretical approach (i.e. linear non-linear dynamics, etc.). We refer the 
interested reader to |28| for a more complete discussion of the implications of 
the existence of a fractal distribution of matter, at least, up to 100h~ l Mpc. 

Let us now remark the predictions for future galaxy redshift surveys. Ac- 
cording to the standard interpretation, the length r ~ 5h~ l Mpc characterizes 
the physical properties of galaxy distributions. Therefore deeper samples like 
CfA2 and SLOAN should simply reduce the error bar, which is now consid- 
ered to be about 10%. A possible variation of r with absolute magnitude, 
due to a luminosity bias, is considered plausible but it has never been quan- 
tified. This should be checked by varying independently absolute magnitude 
and depth of the volume limited samples. However, from this interpretation, 
the value of r = 5h^ x Mpc, corresponding to a volume limited of CfAl with 
M = —19.5, should not change when considering in CfA2 and SLOAN volume 
limited samples with the same solid angle Q and the same absolute magnitude 
limit (M = -19.5). 

In our interpretation, instead, r is spurious, and it scales linearly with the 
radius R s of the largest sphere fully contained in the volume limited samples. 
Therefore we predict for the volume limited sample of CfA2 with M = —19.5 
(with a solid angle of Q = 1.1 sr |16|]) ro ~ lh~ x Mpc (if, in the final version 
of the survey the solid angle is Q = 1.8, the value of R s increases accordingly, 
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Table 3: The volume limited samples of various catalogs (not still published 
and analyzed) are characterized by the following parameters: - Ryi,ihr x Mpc) 
is the depth of the VL sample considered with absolute magnitude limit M VL 
- Q is the solid angle - R s (h~ l Mpc) is the radius of the largest sphere that can 
be contained in the catalog volume. This gives the limit of statistical validity 
of the sample. - ro(h~ 1 Mpc) is the length at which £(r) = 1. (distance are 
expressed in h~ x Mpc). 



Sample 


n (sr) 


RvL 


M VL 


R s 


r 


CfA2 


1.83 


101 


-19.5 


22 


7 


CfA2 


1.83 


160 


-20.5 


36 


12 


SLOAN 


7T 


400 


-19 


185 


60 


SLOAN 


71 


600 


-20 


275 


90 


2dF (South) 


0.28 


550 


-19 


50 


15 


2dF (South) 


0.28 


870 


-20 


100 


30 



and the value of r is shifted up to ~ 9/i _1 Mpc). Note however that for the 
deepest volume limited CfA2 sample (M < —20) we predict instead r$ « 15 4- 
20h~ 1 Mpc. For the volume limited sample of the full SLOAN with M = —19.5 
(Q — 7r), our prediction is that ro ~ 65h~ 1 Mpc. It is clear that however, the 
first SLOAN slice gives smaller values because the solid angle is be small. In 
Tab.|] we report the predictions for r in the next future surveys. 

Finally we would like remark that if it is true that the fractal analysis, 
as good wine, must be taken with "moderation", it is also true that up to 
now the usual approach in this field has suffered of abstinence. We hope that 
the present discussion will evolve in a more constructive and cooperative way 
between the field of statistical mechanics and large scale structure astrophysics, 
hoping to drink together some good wine. 
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